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Abstract 

Germanium is an extremely important mate¬ 
rial used for numerous functional applications 
in many fields of nanotechnology. In this 
paper, we study the crystallization of amor¬ 
phous Ge using atomistic simulations of crit¬ 
ical nano-metric nuclei at high temperatures. 
We find that crystallization occurs by the re¬ 
current transfer of atoms via a diffusive pro¬ 
cess from the amorphous phase into suitably- 
oriented crystalline layers. We accompany our 
simulations with a comprehensive thermody¬ 
namic and kinetic analysis of the growth pro¬ 
cess, which explains the energy balance and 
the interfacial growth velocities governing grain 
growth. For the (111) crystallographic orien¬ 
tation, we find a degenerate atomic rearrange¬ 
ment process, with two zero-energy modes cor¬ 
responding to a perfect crystalline structure 
and the formation of a E3 twin boundary. Con¬ 
tinued growth in this direction results in the 
development a twin network, in contrast with 
all other growth orientations, where the crystal 
grows defect-free. This particular mechanism 
of crystallization from amorphous phases is also 
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observed during solid-phase epitaxial growth of 
(111) semiconductor crystals, where growth is 
restrained to one dimension. We calculate the 
equivalent X-ray diffraction pattern of the ob¬ 
tained nanotwin networks, providing grounds 
for experimental validation. 
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Introduction 

Growth of semiconductor crystals from glassy 
or vapor phases is an extremely important pro¬ 
cess for many applications in (nano) technol- 
ogyjlEl general, crystallization from a dis¬ 
ordered structure is ultimately a diffusive pro- 
cesJ^ and -as such- strongly temperature de¬ 
pendent. Growth, however, is highly suscepti¬ 
ble to the formation of crystal defects, which 
can be copious, mediated by imperfections, 
both related to the environmental variables of 
the physico-chemical treatment ,1^ and to struc¬ 
tural heterogeneities associated with the sub¬ 
strate, such as impurities, lattice mismatch, 
etc.^ Defect formation is typically also temper¬ 
ature dependent and thus a compromise must 
be found to balance reasonable growth rates 
while keeping acceptably-low defect concentra- 
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tions. Since the early times of solid-phase 
crystallization, great emphasis has been placed 
on suppressing the emergence of these imper¬ 
fection^^, as many properties of crystallized 
systems strongly depend on achieving pristine 
structures and a defect-free finish.l^^ 

Among the panoply of possible defects found 
during crystal growth, growth twins stand out 
as one of the more prolihc ones due to low for¬ 
mation energies and a variety of possible gen¬ 
esis pathways In this work we are con¬ 

cerned with twinning in diamond-cubic struc¬ 
tures, such as Si and Ge, whose growth, solidi- 
hcation, and recrystallization have been studied 
extensivelyJ ^ * ^ * ^^ * ^ It is believed that twinning 
may appear as a consequence of internal trans¬ 
formations to rotate the crystal structure to¬ 
wards energetically favorable interfacial orien¬ 
tations However, twinning is also observed 
during epitaxial recrystallization in Si, Ge, and 
their alloys, particularly along the (111) growth 
direction.^ The chief difficulty in studying twin 
nucleation in these systems is that it is gener¬ 
ally impossible to see them form in situ. Sev¬ 
eral authors have hypothesized twins form at 
a very early stage, when the crystal nucleus is 
extremely small. However, twinned overgrowth 
may also be observed forming on a nucleus of 
relatively-large size.^11^ In both of these cases, 
atomistic simulation suggests itself as the ideal 
avenue to unravel the nature of twin nucleation 
and growth. 

Here we report molecular dynamics simula¬ 
tions of Ge recrystallization from amorphous 
substructures at high homologous temperatures 
(fraction of the melting temperature T^). Our 
simulations are framed within the study of 
amorphization-crystallization (a— t-c) processes 
in laser spot heating of GST-basecQ phase¬ 
changing materials (PGMs).^i^^ In GST ma¬ 
terials, the main role of Ge is to accelerate 
the recrystallization process, and thus here we 
study pure Ge as the point of reference for Ge- 
based PCMs.^We start from a critical nucleus 
of crystalline Ge (c-Ge) embedded in an amor¬ 
phous (a-Ge) medium at an initial temperature 
of To = 1100 K, which is representative of the 
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conditions found within the laser spot. We find 
that grain growth is characterized by the for¬ 
mation of intricate twin networks facilitated by 
near zero surface and stacking fault energies 
We accompany our simulations with a full ther¬ 
modynamic analysis to explain the mechanisms 
behind the observed behavior. 

Simulation details 

Molecular dynamics simulations 

We use a Stillinger-Weber potential parameter¬ 
ized by Posselt and Gabriel for Ge,l^ which 
reproduces the experimental values for the 
cohesive energy and lattice constant for the 
diamond-cubic structure, and yields reasonable 
values for the energetics of other crystalline 
phases and the structure of the liquid. Our sim¬ 
ulations are run in the isobaric-isenthalpic en¬ 
semble NpH -where N is the number of parti¬ 
cles, p is the pressure and H is the enthalpy¬ 
using periodic boundary conditions in three di¬ 
mensions. The NpH ensemble was chosen to 
allow for local temperature increases due to la¬ 
tent heat release during the crystallization pro¬ 
cess. To simulate the effect of laser spot heat¬ 
ing, the system is equilibrated to an initial tem¬ 
perature of To = 1100 K, which is approxi¬ 
mately the temperature at the center of the 
spot.^ The starting configuration is generated 
from a perfect diamond structure with lattice 
constant ao = 0.5654 nm, corresponding to the 
value at 1100 K, oriented along the [100], [011], 
and [Oil] directions. The simulation comprises 
62oi X 4402 X 4403 supercell (~35 nm per side, 
where oi = oq, 02 = 03 = oo\/2) containing 
1,920,512 atoms. The central spherical region 
of radius of 2 nm (the nucleus) is then kept 
frozen, while the outer region is melted by fix¬ 
ing the temperature at 3000 K using a Langevin 
thermostat during 100 ps and then quenched. 
At the same time the NpH ensemble maintains 
zero pressure globally. Finally, the entire sys¬ 
tem (the o-Ge block containing the c-Ge nu¬ 
cleus) is further equilibrated at 1100 K during 
an additional 10 ps. The thermostat is then 
turned off during the subsequent crystallization 
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simulations. We have shown that this results 
in a glass transition temperature of approxi¬ 
mately 810 K, in good agreement with labora¬ 
tory experiments for Ge-^^H The procedure just 
described was used by Reina et al.^^ to generate 
planar a/c interfaces to calculate free energies 
and interface mobilities. 

Calculation of the critical nucleus 
size 

In classical nucleation theory, the critical nu¬ 
cleus size r* is governed by the balance between 
the volumetric and interfacial driving forces ex¬ 
pressed, respectively, as the derivative of the net 
free energy release AGq, a^-c and a surface en¬ 
ergy penalty AGg with respect to the radius of 
the nucleus. In principle, this balance must also 
account for the expansion of Ge upon crystal¬ 
lization, which is approximately 10 % less dense 
than its amorphous counterpart (cf. Figure 3b 
in Reina et al.^S) at p = 0 and Tq = 1100 K. 
However, the procedure detailed in the previous 
section to seed an amorphous matrix with crys¬ 
talline grains removes any differential strains by 
construction. This allows us to write: 


crystalline phase at the same temperature is 
Pc = 4.36 X 10^® m“^ (after Reina et al.^. 
From this, Apo.a-i^c ~ —5.07 x 10® J-m“®. 
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Figure 1: Temperature dependence of Apo,a^-c 
expressed on a per atom basis (referred to the 
left vertical axis). Also shown are the internal 
energies u of both amorphous and crystalline 
Ge as a function of temperature (right verti¬ 
cal axis). The heat capacity Cp is calculated 
from the slope of u{T), which for c-Ge results 
in a value of ~ 2.69 x 10“^ eV-K“^ per atom 
(~ 0.36 J-g“bK“^), in excellent agreement with 
experimental measurements.^ 
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AG — AGo,a-i>c + AGs — —^—Apo,a^-c + 47rr^7, 

( 1 ) 

where Apo,a^-c is the volumetric free energy 
density at zero pressure, and 7 is the surface 
(free) energy density, which is orientation de¬ 
pendent: 7 = 7 (d), with 9 representing the 
surface normal with respect to the crystal ori¬ 
entation. We have calculated the atomic free 
energy densities po,a and po,c of the amorphous 
and crystalline phases using thermodynamic in¬ 
tegration, see Reina et dl^ for details. The 
variation of Apo,a^c = go,a - go,c with tem¬ 
perature in units of energy per atom is pro¬ 
vided in Figure When this difference is zero, 
there is phase coexistence, which by definition 
occurs at the melting point, here Tm = 1350 
K. The driving force per unit volume for the 
o —>-0 transformation at zero pressure is read¬ 
ily obtained as: A^fo,a^c = PcAgo^a^c- From 
the figure, at Tq = 1100 K, AgQ^a^c = —0.073 
eV per atom, while the atomic density of the 


As mentioned earlier, AG* is orientation de¬ 
pendent. However, crystalline Ge displays cu¬ 
bic symmetry, which allows us to reduce the ori¬ 
entation space to that contained in the standard 
stereographic triangle whose vertices in the first 
octant are the intersects of the unit sphere with 
the [001], [110], and [111] directions .1^21 Thus, 
we restrict our study of the orientation depen¬ 
dence of 7 ( 6 *) to those three orientation^ The 
inter facial free energies are shown in Figure 
as a function of temperature, where a sur¬ 
face orientation anisotropy can be clearly dis¬ 
tinguished at low temperatures. At Tq however, 
this anisotropy is smeared out by the high ther¬ 
mal diffusivity of the amorphous phase above 
the glass transition temperature, cf. Figure 2 
in Reina et and we find an orientation- 

independent value of 7 ~ 0.08 J-m“^. 

The critical grain size is found by minimizing 

Suitable interpolation schemes within the standard 
triangle can be adopted for a general orientation 9. 
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Figure 2 : Surface free energy as a function of 
temperature for the three surface normals rep¬ 
resenting the vertices of the standard triangle. 


eq. Q: dAG/dr = 0 , which results in 


r 


* 


27 

AQo^q—^c 


( 2 ) 


Replacing 7 and Ago^a^c for their respective 
values, we obtain that r* ~ 0.32 nm. This 
value is approximately 55% of the magnitude 
of the lattice constant Uq and suggests stable 
crystalline grains with only a handful of atoms 
in them, probably implying that very localized 
fluctuations suffice to produce crystal growth 
seeds. 


Interface mobility 

The growth rate of the crystalline phase at the 
expense of the amorphous phase is governed by 
the schematic energy landscape shown in Figure 
IHl The excess atomic flux a^c relative to c— 7 -a 
transitions is governed by a^-c (shown in the 
figure) and results in net interface velocity v{T) 
and grain growth. Mathematically, this can be 
expressed to first order 

v{T) = no exp {-/3Eb) (1 - exp (-/5A^o,a^c)) 

( 3 ) 

where vq is a prefactor and Eb is an activa¬ 
tion energy for the transformation (shown in 
the figure). Eb represents the energy for the 
detachment/reattachment process, which is dif¬ 
fusive in nature. We have devised a special 
procedure to calculate Eb, for which a value 


of 0.42 eV was obtained for the [100] orienta- 
tion^^ at 0 K. A low value of Eb may result 
in faster growth speeds at low temperatures, 
but it also results in faster detachment (a<—c) 
at higher temperatures, where the difference of 
free energies decreases, resulting in lower ef¬ 
fective growth speeds. The interface velocities 



Figure 3: (Free) energy landscape governing 
the crystallization process. 


can also be obtained by direct atomistic sim¬ 
ulation as described by Reina et al.^^. Figure 
m shows results for four distinct orientations at 
To with the (111) being the slowest one, which 
ultimately controls grain growth. 



Angle [°] 


Figure 4: Interface velocity as a function of 
surface orientation at the initial (1100 K) and 
final (1250 K) simulation temperatures. 


Simulation results 

Next we present the simulations of supercriti¬ 
cal grain growth at Tq. While the critical radius 
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calculated earlier suggests a very small stability 
threshold, we have found that -in practice- a 
minimum radius of 2 nm was needed to have 
positive grain growth on the timescales cap¬ 
tured in the MD simulations. The discrepancy 
can be attributed to a number of factors, chief 
of which is the magnitude and frequency of 
thermal fluctuations at these high temperatures 
and small volumes, which lead to low signal-to- 
noise ratios in terms of the stable critical size. 
Other factors such as nonsphericity, and finite 
size effects, may also play a non-negligible role. 
Consequently, in the following we show results 
of 2-nm radius supercritical c-Ge nuclei in an 
amorphous medium. 

Three-dimensional growth of crit¬ 
ical grains 

In a 3D (spherical) nucleus, in principle all 
growth orientations are sampled, which means 
that in materials with sizable interface energy 
anisotropies and/or interface velocities, some 
growth directions will be preferred over others. 
A complete animation of the grain growth pro¬ 
cess starting from (super) critical nuclei is pro¬ 
vided in the Supporting Information. A snap¬ 
shot of the simulation at the point of maximum 
growth -which occurs 9.1 ns after the system is 
equilibrated at Tq- is shown in Figure Atoms 
in the image are colored using structure analy¬ 
sis as implemented in the OVITO visualization 
package,l2Sl which assigns dark blue to atoms 
with diamond cubic structure and orange to 
atoms with hexagonal diamond crystal struc¬ 
ture. 

Closer examination of the atoms with hexag¬ 
onal diamond structure reveals that they be¬ 
long to (111) twin (E3) boundaries, typical of 
the diamond cubic lattice structure. Figure 
shows a region around one such boundary in 
local detail, where the mirror symmetry char¬ 
acteristic of twin plates can be clearly identi¬ 
fied. Moreover, using boundary analysis avail¬ 
able in OVITO, we have examined the atom¬ 
istic structures in Figures |5(a)| and and find 
the emergence of a network of twinned regions, 
as showcased in Figure |5(b)[ The figure shows 
a through-thickness view of the entire grain at 



(b) 


Figure 5: Structure of the crystalline grain 
after 9.1 ns of growth in the NpH ensemble, 
(a) Atomistic structure: Only atoms possess¬ 
ing an ordered crystallographic structure are 
shown. Dark blue and orange spheres represent 
atoms with diamond cubic and twin plate struc¬ 
ture, respectively. Light blue spheres represent 
atoms belonging to the amorphous/crystalline 
interface and dislocation cores. The circular re¬ 
gion in the center of the image indicates the 
extent of the critical grain at the beginning of 
the simulation, (b) Through-thickness view of 
the nanotwin structure corresponding to (a). 

the exact same time as Fig. |5(a)[ As depicted, 
in the lower half of the grain twins are elongated 
along (211) directions, while in the upper half 
a three-dimensional arrangement is formed. 

As discussed earlier, grain growth is con¬ 
trolled by the magnitude of the driving force 
and the interface mobility. Both of these quan¬ 
tities are temperature dependent. The exother¬ 
mic nature of the a —>-0 reaction (i.e. internal en¬ 
ergy density difference Aua^c < 0) results in a 
local energy deposition that increases the global 
system temperature. Thus, the temperature of 
the system correlates directly with the volume 
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Figure 6: [Oil] view of the atomistic configu¬ 
ration of a twinned region. 


of material transformed. Such correlation is 
clearly visible in Figure where both quanti¬ 
ties are exactly proportional to each other with 
a proportionality constant of ~ 8.01 x 10“^^ 
m^ K“h The temperature is seen to increase 
from To to a final value of approximately 1250 
K. This effectively arrests the growth process, 
as dictated by the sharp decrease in mobility 
at such temperature (cf. Figure]^. We remark 
that this arrest is partially an artifact of the 
simulations which limits unrestricted heat flow 
due to periodic boundary effects. 

Analysis of growth of a/c bicrys¬ 
tals 

The appearance of twins during the growth 
stage of critical Ge grains may obey the energet¬ 
ics of two different scenarios. On the one hand, 
there is ample evidence in the literature that 
in Ge twinning emerges if the growth direction 
deviates appreciably from the preferred growth 
direction, understood as that which results in 
the lowest interfacial energy. Twins would then 
appear to alter the internal crystal orientation 
and bring it closer to the preferred one.^^^ How¬ 
ever, other works have pointed out that if a 
new facial orientation was the only advantage 
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Figure 7: Evolution of the temperature with 
time during growth of a critical nucleus under 
Nph simulations conditions. The crystallized 
volume is also shown as a function of time. Se¬ 
lected snapshots of the grain structure corre¬ 
sponding to three distinct instants (2, 5, and 8 
ns) are shown. 


gained by twinning, twinned crystals should not 
be much larger than twice the size of a single 
crystal, something at odds with observations of 
twinned crystals being more than ten times as 
large as untwinned ones.^ This hypothesis is 
strongly weakened by the lack of a noticeable 
surface energy anisotropy at 1100 K according 
to our calculations (cf. Figure [^. 

The alternative scenario is that twins are 
a manifestation of a growth mode that relies 
on the indistinct formation of ordered atomic 
layers with the correct stacking sequence and 
stacking faults. This is the same growth mode 
observed under the so-called solid-phase epitax¬ 
ial recrystallization (SPER) process of {111}- 
oriented crystals, as well as by liquid epitax¬ 
ial growth of crystals with the same orientation 
via chemo-physical vapor deposition. There is 
ample evidence of twin formation in the litera¬ 
ture for both of these processes in Ge, partic¬ 
ularly at high temperatures .1231221 This mecha¬ 
nism is controlled by the value of the stacking 
fault energy 7 sf| 3 which ranges between 0.07 

^Twin boundaries correspond to ‘half’ a stacking 
fault and so the twinning propensity correlates directly 
with the stacking fault energy 
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and 0.09 J-m“^ according to several measure¬ 
ments.®®^ By contrast, the interatomic po¬ 
tential employed in our simulations predicts 
zero stacking fault energy.®^ Evidently then, the 
model for Ge employed here offers no imped¬ 
iment to the favorable formation of epitaxial 
twins. However, while at low temperatures this 
might clearly result in an overestimation of the 
volume fraction of stacking faults and/or twins 
when conditions are conducive to their forma- 
tion,®®! it is reasonable to assume that val¬ 
ues of 7 sf on the order of the experimentally- 
measured ones result in zero effective stacking 
fault energy at a temperature of 1100 K via 
thermal softening. 

To ascertain which mechanism is responsible 
for the observed formation of nanotwin net¬ 
works, next we carry out MD simulations of 
o/c bi-crystals at Tq oriented along three se¬ 
lected directions: [111] (low mobility, cf. Fig¬ 
ures and 1^, and [100] and [110] (high mo¬ 
bility). These are qualitatively similar to other 
simulations of the SPER process using atom¬ 
istic methods.®®^ The three surface orienta¬ 
tions simulated here are schematically shown 
in Figure relative to a [110] view of the Ge 
diamond cubic lattice. 



Figure 8: Schematic view along the [lIO] direc¬ 
tion of the Ge diamond cubic lattice, with the 
(111) plane highlighted. Image obtained with 
Wolfram CDF Player.®^ 


[Ill] amorphous/crystalline bi-crystals 

The starting microstructures (o/c bi-crystals) 
are generated in the manner described by Reina 
et al.® The computational cell has dimen¬ 
sions of 20.8 X 19.4 X 39.2 nm containing 698,880 
atoms. The system is again equilibrated at 1100 
K and let to evolve in the NpH ensemble. Two 
animations illustrating the process are provided 
in the Supporting Information. Growth of the 
crystalline phase proceeds via the formation of 
an intricate twin network, an image of which is 
shown in Figure [^ Twin boundaries are shown 
as green-colored surfaces, and are seen to form 
a quasi-hexagonal network as demanded by the 
topological structure of a set of interconnected 
{111} surfaces. The structures are reminiscent 
of coral-like porous networks in synthesized ce¬ 
ramic materials 



Figure 9: Twinning network formed from the 
growth of a crystalline Ge half-crystal along the 
[111] direction. 


[100] and [110] amorphous/crystalline bi¬ 
crystals 

The size of the computational cells employed 
to study growth along the [100] and [110] direc¬ 
tions was 20.4 X 20.4 x 39.6 nm with 725,760 
atoms, and 20.4 x 20.8 x 40.0 nm with 748,800 
atoms, respectively. As shown in the corre¬ 
sponding animations (Supporting Information), 
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crystallization along these directions results in 
growth of a homogeneous Ge crystal, forming 
essentially no defects. Here, atom rearrange¬ 
ments from the amorphous into the crystalline 
phase occurs by forming atomic planes with the 
correct stacking sequence. 


XRD analysis 


An important part of the analysis of the simula¬ 
tions is their experimental validation. Ge crys¬ 
tals can be examined by a variety of sources, 
from transmission electron microscopy (TEM), 
to Raman spectroscopy (RS), and X-ray diffrac¬ 
tion analysis (XRD). However, twin boundaries 
are not sources of strain and are thus difficult 
to detect via conventional TEM analysis. In 
contrast, they act as scattering agents to X- 
rays and do leave an imprint on diffraction pat¬ 
terns. Therefore, we have calculated the equiv¬ 
alent XRD signature for pure crystalline and 
amorphous samples, as well as for the multi- 
twinned structure shown in Eig. using the 
Debyer code®' considering a X-ray source with 
wave length of 1.542 A at 0 K. The resulting 


pattern is shown in Eigure where the in¬ 
tensity peaks represent the different scattering 
directions. The figure reveals clear differences 
in the footprints of the three structures consid¬ 
ered, namely, no structure for the amorphous 
system, well marked peaks for the ideal crys¬ 
tal, and softened peaks for the twinned crys¬ 
tal. The XRD pattern showed in the Figure 
is in excellent agreement with experimental re¬ 
sults for pure crystalline Ge but only in mod¬ 
est agreement for amorphous Ge.®' This may 
indicate that the generated amorphous struc¬ 
tures may not be fully optimized in terms of 
their atomic configuration, likely a result of us¬ 
ing unphysically-high heating and cooling rates 
to entrap a liquid structure into a disordered 
solid. With regard to multitwinned structures. 
Fig. provides a pathway for their detection 
in future experiments. 


Discussion and conclusions 

Although it is clear from the literature that 
twinned Ge crystals may emerge during crys¬ 
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Figure 10: Diffraction patterns for amorphous, 
perfect crystal and twinned Germanium. The 
atomic structures were minimized before be¬ 
ing analized with Debyei®' considering a X-ray 
source with wave length of 1.542 A. 


tallization at low temperatures to favor low en¬ 
ergy interface orientations (and thus decrease 
the critical nucleus size), our simulations con¬ 
clusively show that the origin of the twinning 
network observed during crystallization of Ge 
grains from amorphous structures at Tq lies in 
the energetic degeneracy observed for the stack¬ 
ing sequence of {111} planes. Although this 
effect is favored by construction in our simula¬ 
tions (due to a zero stacking fault energy pre¬ 
dicted by our atomic model), the overall effect 
of 'jsF in materials such as Ge at these high 
temperatures is likely to be negligible in any 
case. The result is the spontaneous formations 
of multiply twinned structures along each of the 
three equivalent (111) directions. 

The other notable observation is that grain 
growth at 1100 K is controlled by low inter¬ 
face mobilities. This together with a small 
critical radius for stable crystalline nuclei, sug¬ 
gests the development of nano crystalline or 
very fine-grained structures, as is indeed the 
case experimentally. The reason for this is 
that 1100 K is near the tipping point where 
the mobility sharply decreases from its max¬ 
imum value. This is compounded by latent 
heat deposition released during the exothermic 
a^c process, which increases the temperature 
beyond that tipping point. The phenomenon 




















where crystallization is fueled by the intrinsic 
latent heat release is well known and referred to 

as explosive crystallization.l^^lil Pqj- this, how¬ 
ever, an increase in temperature should re¬ 
sult in growth acceleration by a surge in in¬ 
terface velocity, which is not the case in the 
temperature regime considered here. In all, 
high nucleation rates due to small critical radii, 
slow mobilities due to high temperatures, plus 
high-twinning propensities result in the no¬ 
toriously fine-grained nanostructures reported 
for laser-induced Ge crystallizationWe 
have recently proposed a thermodynamically- 
consistent phase field model to predict these 
microstructures .1221 However, intrinsic twinning 
was not a feature of those simulations and we 
believe that the present atomistic simulations 
provide a new piece of physics that must be in¬ 
corporated into such higher-level models. 

Twin boundaries may also act as scattering 
agents for elastic and electromagnetic waves, 
and may impact the value of fundamental con¬ 
stants such as the thermal conductivity or elec¬ 
tric susceptibility. Indeed, it has been observed 
that the appearance of twins during epitaxial 
growth of Si wafers resulted in faulty devices, 
while for other defects, such as e.g. extrinsic 
stacking faults, it did not®® (albeit perfect 
Ge crystals have been grown in the [ 111 ] di¬ 
rection as well®). This may be of importance 
in GST materials where high contrast between 
amorphous and crystalline phases in terms of 
these properties is desired. 
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